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Abstract. Yagi-Uda antennas are known to be difficult to design and 
optimize due to their sensitivity at high gain, and the inclusion of nu- 
merous parasitic elements. We present a genetic algorithm- based auto- 
mated antenna optimization system that uses a fixed Yagi-Uda topology 
aind a byte-encoded antenna representation. The fitness calculation al- 
lows the implicit relationship between power gain and sideiobe/backlobe 
loss to emerge naturally, a technique that is less complex than previous 
approaches. The genetic operators used are also simpler. Our result- 
s include Yagi-Uda antennas that have excellent bandwidth and gain 
properties with very good impedance characteristics. Results exceeded 
previous Yagi-Uda antennas produced via evolutionary algorithms by at 
least 7.8% in mainlobe gain. We also present encouraging preliminary 
results where a coevoiut ionary genetic algorithm is used. 


1 Introduction 

Automated antenna synthesis via evolutionary design has recently garnered 
much attention in the research literature [12]. Underlying this enthusiasm is an 
issue that many designers reaulily acknowledge - good antenna design requires 
not only knowledge aind intelligence, but experience amd artistry. Thus automat- 
ed design techniques and tools have been lauddng. Evolutionary algorithms show 
promise because, among search algorithms, they au:e able to effectively search 
large, unknown design spaces. 

The particular antenna we study in this paper is the Yagi-Uda, first proposed 
in 1926 [14]. We chose this type of antenna because it presents diflScult design 
amd optimization chadlenges, and because it was previously studied with respect 
to. evolutionary design [7]. The Yagi-Uda antenna is comprised of a set of paLraJIel 
elements with one reflector element, one driven element (driven from its center), 
amd one or more director elements (see Fig. 1). The highest gaun can be achieved 
ailong the axis and on the side with the directors. The reflector element reflects 
power forwards and thus acts like a smadl ground plane. The design pauauneters 
consist of element lengths, inter-element spacings, and element diameters. 
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Fig. 1. Typical Yaigi-Uda antenna. 


The application that we use is taken from [7). It involves designing a spe- 
cial feed for the Arecibo 305- meter spherical reflector in Puerto Rico [3]. The 
antenna was to be used to search for primeval hydrogen having a redshift of 
approximately 5. Neutral hydrogen line emission is at a frequency of 1420 MHz; 
thus the frequency region of interest was about 235 MHz, Preliminary studies 
indicated that the band from 219 to 251 MHz was of the greatest interest, par- 
ticularly from 223 to 243 MHz. The most important design goal was for the 
feed to have sidelobes/backlobes at least 25 dB down from the mainbeaun gain 
in the region from 70® < <p < 290°, due to the interference which came from 
surrounding radio and TV towers. Of lesser importance was that the E-plane 
(the plame paradlel to the plaine of the antenna) amd H- plane (perpendicular to 
the E-plane) beamwidths be about 50°. 

Voltaige Standing Wave Ratio, or VSWR, is a way to quantify reflected- 
wave interference, and thus the amount of impedance mismatch at the junction. 
VSWR is the ratio between the highest voltage and the lowest voltage in the 
signal envelope along a tramsmission line [13]. The VSWR was desired to be less 
than 3 and the gaun was to be maodmized, limited by the wide beaunwidth. The 
feed would be mounted over a 1.17 meter square ground plane-that is, a ground 
plane only 0.92 A in size. 


2 Antenna Representation and Operators 

The representational scheme used is similar* to that taken from [7]. As shown in 
Fig. 2, this scheme is comprised of 14 elements, each one encoding a length and 
speicing value. Eatch floating point vadue was encoded as three bytes, yielding 
a resolution of 1/2^"* per value. The first pair of values encoded the reflector 
element, the second paur encoded the driven element, and the remaining 12 pairs 
encoded the directors. One point crossover was used with cut points allowed 
between bytes. Mutation was applied on individual bytes. 
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Fig. 2. Genetic representation of a 14-eIement Yagi-Uda antenna. 


Fladius values were constrained to 2, 3, 4, 5, or 6 mm. All elements within 
a given individual were assigned the same radius value. Element lengths were 
constrained to be symmetric around the j-axis and between 0 and 1.5 A. Ele- 
ments having zero length were removed from the antenna; as a consequence, a 
constructed antenna could have less than 14 elements. Spacing between adjacent 
elements (aJong the z 2 Lxis) was constrained to be between 0.05 A amd 0.75A. The 
wavelength A was 1.195 meters, the wavelength of 235 MHz. 

3 Experimental Setup 

Experiments were set up as follows. The NEC2 simulation program [4] was used 
to evaluate all antenna designs. We used a parallel master/slave generational 
genetic adgorithm with a population size of 6000. One point crossover across 
byte boundaries was used at a rate of 80%. Mutation was uniform across bytes 
at a rate of 1%. Runs were executed on a 32-node Beowulf computing cluster [111. 

The wire geometry encoded by each individual chromosome was first translat- 
ed into a NEC input deck, which was subsequently sent to the NEC2 simulator. 
The segment size for all elements was fixed at 0,1 A, where A was the wavelength 
corresponding to 235 MHz. The source element for excitation was specified to be 
the middle segment of the driven element. The z location of the reflector element 
was always set to 0. The antenna was analyzed in free space. 

The simulator was instructed to sample the radiation pattern of each indi- 
vidual at three different frequency values: 219, 235, and 251 MHz, representing 
a 13.6% bandwidth. Batch radiation pattern was cadculated at set to 0® aind 6 
varying between 0® and 355®, the latter sampled at 5® increments. VSWR values 
were adso calculated for each of the three frequencies. 

Fitness was expressed as a cost function to be minimized. The calculation 
was as follows: 


where: = lowest gaun of all frequencies measured at ^ = 0® atnd 0 = 0®, K = 

VSWR at the ith frequency, and 

^ / 0.1 if K < 3 

^ \ 1 if K > 3 





Lackiru; from Miis w;ls i tv.rm involving sldrlofx'/hai klobe attenua- 

tion. VVf? chose not iih lude such a term l)erause wr* re.isoned tliat i\a the mainlobe 
j^ain mcre<used, tlie si<l<*lohes/l);u klobes wr>uid decroiLso in size. 


4 Experimental Results 


Thirteen runs were executed untier iliffering random number streams for com- 
parison purposes. Table I summarizes the run data for the best antenna found 
in each run of 100 generations. Fig. 3 shows the radiation pattern from the best 
antenna found (run 13). It exhibits 10.58 dB and hcis a VSWR of 2.02 at its 
center frequency. Its sidelobe/backlobe gain at this frequency is 3.07 dB. Fig. 4 
shows a diagram of the antenna's physical structure. 

To increase simulation speed, the evolved antennas were produced without 
the presence of a ground plane - an idealized setting. Adding a ground plane 
thus simulates more realistic conditions. We removed the reflector element and 
simulated the best antennas found over a ground plane of L17 meters [7]. We 
found the performance increased - at the center frequency the mainlobe gain 
was 12.52 dB and the VSWR was 2.39. At 291 MHz, the gain was 11.33 dB, and 
at 251 MHz, the gain was 11.15 dB. In contrast, the antenna produced in [7] 
exhibits gains of 10.36, 10.91, 10.34 dB at 219, 235, and 251 dB, respectively. 
Thus the antenna from run 13 has a minimum performance increase of 7.8% as 
compared to the previously reported antenna. 
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Table 1. Results from the best individual after 100 generations for each of the 13 runs 
(dB is measured at = 0°, ^ = 0^). 







































219 MHz 235 MHz 251MHz 



Fig. 3. Radiation pattern of the best evolved antenna without a ground plane, mea- 
sured at 0"^ < 0 < 360°, 0 = 0°, for 219, 235, and 251 MHz, respectively. (The scale is 
2 dB per division. Inner ring is -12 db, outer ring is 12 dB.) 
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Fig. 4. The best Yagi-Uda antenna from run 13. The radius of aJl elements was 3 mm. 


5 Coevolutionary Algorithm - Method and Preliminary 
Results 


A cnevolutionary algorithm wi\s also to th(j iuitenna optimization 

prol)lern doscTihprl abov(\ The experiments are oni^oiiig ;ks of this writing, an<l we 
briefly mention some encouraging initial results. The algorithm used is similar to 
that presented in [9[. Two popiilations are used: one consisting of antenna designs 
.as described above, and one consisting of target vectors. The fundamental idea 
is that the target vectors encapsulate level-of-difficulty. Then, under the control 
of the genetic algorithm, the target vectors evolve from easy to difficult based 
on the level <if proficiency of the antenna population. 

Each target vector consists of a set of objectives that must be met in order 
for a target vector to be ‘‘solved.’’ A tau^get vector consisting of two values: the 
mainlobe gain (in dB) and a VSWR value. A target vector was considered to be 
solved by a given antenna if: 

^target < ^;arget > 

where Gi, is lowest gain of all frequencies measured at ^ = 0® and ^ = 0®, and 
Vl is lowest VSWR of all frequencies. For example, an antenna with a Gi value 
equal to 5 dB and a value equal to 8 would solve the target vector (2, 12) but 
not (7,12). 

V 2 flues for tau*get gaun ramged between 0 dB (eaisy) and 12 dB (difficult). 
Tairget VSWR values ranged between 12 (easy) amd 3 (difficult). Target vectors 
atre represented as a list of floating point vaJues that au’e mutated individually by 
randomly adding or subtracting a small aunount (5% of the largest legal value). 
Single point crossover was used, amd crossover points were chosen between the 
vadues. 

The general form of the fitness caJculations are from [9], In summary, anten- 
nais are rewarded for solving difficult target vectors. The most difficult target 
vector is defined to be the target vector that only one antenna can solve. Such 
a taurget vector gau*ners the highest fitness score. Target vectors that are unsol v- 
abie, or axe very eaisy to solve by the current amtenna population, are given low 
fitness scores. 

We ran our coevolutionairy aJgorithm for 200 generations using 1600 indi- 
viduals in both populations. In the amtenna population, crossover and mutation 
rates were 0.8 amd 0.1, respectively. In the target vector population, crossover 
amd mutation rates were 0.8, 0.5, respectively. 

The highest- fit ness individuad catme from generation 199. It had mainlobe 
gains of 8.30, 8.51, and 8.30 dB at 219, 235, and 251 MHz, respectively. While 
performamce is less tham the runs from above, it was achieved with a much 
smaller population, and it is currently our single data point. 

Fig. 5 shows a plot of how the highest fitness target vectors varied during the 
run. Such plots can give insight regarding the difficulty of achieving one objective 
at the expense of another. In the plot, we see that difficult VSWR leveb (near 
3.0) are attainable eau^ly on and remadn so throughout the run. The algorithm 
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Fig. 5. Coevolution run: plot of gain and VSWR for the best target vector over 200 
generations. 


6 Discussion 

Small improvements in antenna performance can be significant in many ap- 
plications. Because of their numerous design variables, complex behavior, aind 
sensitivity to parauneters, Yetgi-Uda atntennas are notoriously difi&cult to opti- 
mize. Our experiments produced several excellent antennas in a relatively small 
number of generations. When simulated over a finite ground plane, the highest 
performaince antenna found exhibiting a mainlobe gain that was 7.8% higher 
than a previously-reported antenna. 

Previous work has explicitly included a sidelobe/backlobe term in the fitness 
function in order to minimize radiation outside of the desired direction [7]. We 
did not include an explicit sidelobe/backlobe term but rather relied on the fact 
that the radiation pattern of an antenna is a zero sum quantity - increasing the 
intensity in one direction will implicitly reduce the amount of radiation in other 
directions. 

Finally, we are encouraged by our preliminary results produced using coevo- 
lutionary optimization. There we saw an antenna generated that had very good 
properties while requiring less evaluations than the standard GA approach. 
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